Spin-orbit interactions and spin-currents from an exact-exchange Kohn-Sham method 
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An exact-exchange spin-current Kohn-Sham method to treat non-collinear spin, magnetic effects, 
currents, spin-currents, and spin-orbit interactions self-consistently on equal footing is introduced. 
Spin-orbit interactions are shown to induce spin-currents. Results for silicon and germanium are 
presented. 



The electronic structures of crystalline materials con- 
taining heavier elements are strongly affected by spin- 
orbit effects. Even in the band structure of germanium, 
a third row element, spin-orbit splittings of about 0.25 
eV occur. Thus for an accurate description of electronic 
structures of most interesting materials a treatment of 
spin-orbit effects is highly desirable, ff not just sim- 
ple bulk semiconductors but also other materials, e.g., 
oxides or metals or one- and two-dimensional periodic 
systems like molecular wires or surfaces, shall be con- 
sidered then in addition to spin-orbit effects also spin- 
polarization, including non-collinear spin-polarization, 
needs to be treated, ff magnetic properties shall be stud- 
ied, an inclusion of currents is necessary, fn the emerg- 
ing field of spintronics besides density currents also spin 
currents may be of importance. Materials with proper- 
ties of technical interest, e.g. for new semiconductors 
or magnetooptical storage devices or for spintronics, of- 
ten are characterized by an interplay of spin-orbit inter- 
actions, non-collinear spin-polarization, orbital currents, 
and spin-currents. Therefore it seems highly desirable 
to develop an electronic structure method that treats 
all of the mentioned effects in a unified way on equal 
footing. This work outlines such a method within the 
framework of density-functional theory (DFT) and, as a 
central result, shows how spin-orbit interactions gener- 
ate spin-currents and how the latter can be treated in 
a self-consistent manner. As an illustration of the new 
method we calculate spin-orbit effects in the bandstruc- 
tures of silicon and germanium with special emphasis on 
spin-currents induced by spin-orbit effects. 

DFT is applied almost exclusively within the Kohn- 
Sham (KS) formalism^. In recent years KS methods that 
treat the exchange energy as well as the local multiplica- 
tive KS exchange potential exactly were developed 2 ' 3,4,5 . 
Such exact-exchange (EXX) Kohn-Sham (KS) methods 
proved to describe the electronic structure of periodic 
systems, in particular bandstructures and band gaps of 
semiconductors, distinctively better than conventional 
KS methods based on the local density approximation 
(LDA) or generalized gradient approximations (GGAs) 2 . 

EXX KS methods not only yield improved bandstruc- 
tures. Furthermore, they represent a systematic im- 
provement over LDA and GGA methods in the sense 
that the largest fraction of the exchange-correlation func- 



tionals for energy and potential, i.e., the exchange part, 
no longer needs to be approximated but is treated ex- 
actly. This is crucial for the treatment of magnetic effects 
within DFT. Already about 20 years ago Vignale and Ra- 
solt introduced current-spin-density-functional theory^ 
which, in principle, enables a treatment of both spin 
and magnetic effects. In practice, however, current- 
spin-dcnsity-functional theory could not be applied to 
realistic systems because reliable and tractable approxi- 
mate current-density functionals were not available. The 
EXX framework offers the opportunity to avoid at least 
exchange approximations and thus opens the route to 
current density functional methods^. Indeed, we re- 
cently presented an EXX spin-current density-functional 
method that is based on a spin-current density-functional 
theory (SCDFT)2. SCDFT generalizes the current spin- 
density- functional theory of Vignale and Rasolt by taking 
into account spin-currents in addition to the quantities 
considered in current spin-density- functional theory, i.e., 
the electron density, the x-, y-, z-components of the spin- 
density, and the x-, y-, z-components of the paramagnetic 
current density. The spin-currents that are additionally 
considered in SCDFT are the x-, y-, z-components of the 
paramagnetic currents of the x-, y-, z-components of the 
spin-density, i.e., paramagnetic currents of the magneti- 
zation. 

The EXX SCDFT method enables the treatment 
of non-collinear spins, currents of the electron density 
and, going beyond currcnt-spin-density-functional the- 
ory, spin-currents. By introducing in this work spin-orbit 
and scalar relativistic effects in the EXX SCDFT method 
the desired comprehensive description of electronic and 
magnetic properties of materials including relativistic ef- 
fects, in particular spin-orbit effects, becomes possible. 

Our starting point, the EXX SCDFT method of Ref. 
, employs basis sets of plane waves and thus requires to 
treat core electrons via pseudopotentials. Relativistic ef- 
fects including spin-orbit splittings are generated almost 
exclusively in the region close to nuclei. In this region or- 
bitals and subsequently the electronic structure is deter- 
mined by the pseudopotential. In order to introduce rela- 
tivistic effects, we therefore introduce spin-orbit-resolved 
normconserving EXX pseudopotentials vf p with £ de- 
noting the spatial orbital momentum and j = I ± 1/2 
(j = 1/2 for £ = 0) denoting the total, i.e., spatial 
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plus spin, angular momentum. Instead of nonrelativis- 
tic ^-dependent EXX s-,p-,d-pseudopotentials we thus 
construct ^'-dependent EXX Si/2-,,Pi/2- #3/2-^3/25^5/2- 
pseudopotentials that take into account all relativistic 
effects and employ them as separable pseudopotentials 



sep 



(r,r') = 



3 



m=-j 



x ^(^(^(^^(^(f). (1) 

In Eq.JQ) the £lij m denote central field spinors, i.e., two- 
component angular momentum eigenfunctions with to- 
tal angular momentum j and corresponding magnetic 
quantum number m that result from angular momentum 
coupling of spin eigenfunctions with quantum numbers 
s = 1/2 and m s — ±1/2 represented by two dimensional 
unit vectors and spherical harmonics with quantum num- 
bers £ and mi- The pseudopotentials v pp are constructed 
in such a way that the valence orbitals resulting from a 
pseudopotential EXX calculation for the atom, i.e., the 
atomic pseudoorbitals <pi j , are identical to the large com- 
ponents of the valence orbitals from a fully relativistic all 
electron EXX calculation beyond the radius r c and have 
the same eigenvalues. 

The KS Hamiltonian operators of the pseudopotential 
EXX calculation for the atom and of the final plane wave 
EXX-SCDFT calculation are nonrelativistic except that 
they contain the pseudopotentials v pp . The latter in- 
troduce the relativistic effects in particular the spin-orbit 
interaction. As a result the plane wave EXX-SCDFT cal- 
culation then differs from a standard nonrelativistic pseu- 
dopotential EXX calculation not only by the fact that the 
orbitals are two-dimensional spinors that contain spin- 
up and spin-down components which, in general, both 
are nonvanishing, but in addition by the coupling of the 
two spin components through the the £j-dependent pseu- 
dopotentials. The latter are obtained along the lines de- 
scribed in Ref. nam , however, without averaging the 
two Zj-dependent pseudopotentials v pp with j = t— 1/2 
and j =£+ 1/2. 

Technically we construct Zj-dependent pseudopoten- 
tials v pp for a fully relativistic four-component pseu- 
dopotential calculation. The corresponding pseudoor- 
bitals occuring in the separable pseudopotentials v s ^ p 
of Eq. at first, are four-component spinors. 

We then use only the renormalizcd large components 
<t>ij of the original four-component pseudoorbitals to 
construct separable pseudopotentials v^ p for a two- 
component valence-electron calculation that takes into 
account relativistic effects only through the pseudopo- 
tentials. Strictly speaking this is somewhat inconsistent. 
However, the charge of the pseudoatoms is very low, for 
Si and Ge just 4 a.u. and it is save to assume that di- 
rect relativistic effects, i.e., relativistic effects that do not 
have their origin in the pseudopotential, but would result 
from a fully relativistic four-component treatment of the 
valence electrons, are negligible for such small effective 



nuclear charges. Indeed for the atoms Si and Ge the ener- 
getic differences from a two-component pseudopotential 
calculation and a fully relativistic four-component pseu- 
dopotential calculation are of the order of /iHartree and 
thus completely negligible. 

In order to analyze how relativistic effects are taken 
into account in a EXX SCDFT calculation via the j- 
dependent pseudopotentials v pp constructed with the 
relativistic EXX approachiSiii, we consider for the mo- 
ment the KS orbitals of a bulk system as linear combina- 
tions of atomic pseudoorbitals 4>ij . One effect of the rela- 
tivistic pseudopotentials is that s- and p-orbitals are con- 
tracted and lowered in energy while d- and /-orbitals are 
more diffuse and raised in energy. The influence of these 
effects on the KS orbitals in the bulk are automatically 
taken into account in a standard EXX treatment. The 
/-dependence of the pseudopotentials, which is a conse- 
quence of spin-orbit interactions, on the other hand, as 
mentioned above, leads to two-dimensional spinors as or- 
bitals and causes features not present in standard EXX 
methods, like spin-currents. 

For a more detailed discussion of the effect of spin- 
orbit interactions and in particular of the generation of 
spin-currents we consider a two-component all-electron 
Pauli-Hamiltonian operator of an electronic system in a 
magnetic field including spin-orbit interactions: 



H= f + V ee + v ext +H ma s + H so 
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f + V ee + Jdr S T V(r) J(r) . 
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In Eq. 0, T and V ee are the operators of the kinetic 
energy and the electron-electron repulsion, v ext is the 
operator generated by the external electrostatic poten- 
tial v ext , usually the potential of the nuclei, H ma 9 are 
the parts of the Hamiltonian operator caused by a mag- 
netic field B with accompanying vector potential A, and 
H so describes the spin-orbit interaction. By the po- 
sition of the i-th electron is denoted, by Pi the corre- 
sponding canonical momentum operator. The vector a 
contains the Pauli spin matrices. The sum in the first 
line of Eq. [2] runs over all N electrons. The vec- 
tor £ has four components, £0 being a 2x2 unit ma- 
trix and Ei, £2, and £3 being the Pauli spin matrices 
er x , cr y , and cr z . The four components of the vector 

J(r) are the density operator X)i=i — r i) = ^o( r ) 
and the x-, y, z- components of the current operator 
(2) Yh=i Px,i 5(r-ri) + <S(r-r;) p x , t = Ji(r) etc. The 
4x4 matrix V(r) is composed of matrixelements V^„(r) 
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with /i, v — 0, 1, 2, 3 and is given by 

fv ext (r) + ^ A x (r) A y {v) A z (r) \ 
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In Eq. [jj], u e xt,x( r ) = dv ext (r)/dx. The fact that the 
spin-orbit interaction Hamiltonian operator _ff 50 can be 
expressed according to the last line of Eq. [2] together 
with Eq. |3] is crucial, because this represents the ba- 
sis for the SCDFT recently introduced in Ref. 9]. The 
spin-orbit interaction only leads to additional terms in 
the matrix V of Eq. [3] , which do not interfere with the 
derivation of the SCDFT of Ref. 0, i.e., the latter deriva- 
tion also holds in the presence of the spin-orbit Hamilto- 
nian operator H . 

According to Ref. the many-electron KS equation 
associated with a Schrodinger equation with an Hamilto- 
nian operator given by Eq. [2] is determined by the KS 
Hamiltonian operator 

H s = f + Jdr S T V s (r) J(r) (4) 

with the 4x4 matrix V s representing the KS potential, 
which is given by V s (r) = V(r) + U(r) + V xc (r). The 
Hartree potential U and the exchange-correlation poten- 
tial V xc can be represented by 4x4 matrices. Within the 
matrix U only the component Uqq(t) = Jdr' ^_r' |^ > the 
standard Hartree potential, is different from zero. The 
components ~V XC ^ U of the exchange-correlation potential 
V IC are given by 

WW = 7~7T ■ (5) 

<VV( r ) 

Within the considered SCDFT all 16 components of 
the spin-current density p of the KS system equal those 
of the real electronic system. The components p^ v 
with p, v = 0,1,2,3 of p are defined as follows: poo is 
just the regular ground state electron density, p^o with 
fj, = 1,2,3 represent the x-, y-, and z-components of the 
spin-density, pq v with v = 1,2,3 represent the x-, y-, and 
z-components of the paramagnetic current of the electron 
density, while p^ with p,, v = 1,2,3 represent the com- 
ponents of the paramagnetic currents of the spin-density, 
i.e., currents of the magnetization. 

In the absence of spin-orbit interactions a correct de- 
scription of all 16 components of the spin-current den- 
sity p through the KS system is desirable for systems 
that exhibit non-vanishing currents of the spin-density, 
similarly as it is desirable in standard, i.e., non-current, 
density-functional theory to treat spin-polarized systems 
within spin-density-functional theory in order to cor- 
rectly describe both spin-up and spin-down electron den- 
sities. However, in the absence of spin-orbit interac- 
tions, the Hamiltonian operator Eq. @ does not cou- 
ple to currents of the spin-density, i.e., does not couple 



to p^u with p,v = 1,2,3. Therefore it is possible in 
these cases to apply the standard current spin-density 
formalism of Vignale and Rasolt, which considers only 
the ground state electron density, pod the spin-density 
p M o with p = 1,2,3, and the paramagnetic current of 
the electron density p$ v with v = 1,2,3, similarly as it 
is possible to treat spin-polarized systems via non-spin- 
polarized density-functional theory although this might 
affect the accuracy in practical application, which re- 
quire the use of approximate density- functionals. In the 
presence of spin-orbit interactions the situation is differ- 
ent. Because the Hamiltonian operator [2] now couples 
to components of the paramagnetic currents of the spin- 
density it is mandatory that the KS system correctly de- 
scribes all 16 components of the spin-current density p MJ/ , 
i.e., the standard current spin-density functional formal- 
ism is no longer applicable in the presence of spin-orbit 
interactions. 

In our case we do not treat spin-orbit interactions by 
the explicit spin-orbit term H so of Eq. (J2J, but via the 
j-dependent pseudopotentials v s ^ p of Eq. (Q. However, 
j-dependent pseudopotentials, by construction, generate 
an effect corresponding to that of the spin-orbit term 
H so and thus like the latter generate spin-currents and 
couple to them. It is important to notice that spin- 
currents, due to the continuity of orbitals, are present 
not only close to the nuclei, but also in the interstitial 
region between the atoms despite the fact that they are 
generated close to the nuclei by pseudopotentials. Thus 
spin-currents affect the bonding and generally the elec- 
tronic structure. A full and self-consistent treatment of 
the effects of spin-currents in DFT therefore requires a 
spin-current density-functional treatment as provided by 
the EXX SCDFT method of this work. 

Note that an explicit treatment of spin-orbit inter- 
actions via H so in a selfconsistent KS method would 
be problematic, because the spin-orbit term H so aris- 
ing from the relativistic many-electron Hamiltonian op- 
erator via the singular Wouldy-Wouthausen transforma- 
tion leads to severe instabilities. Two-component meth- 
ods that do not employ pseudopotentials therefore treat 
spin-orbit effects either perturbativelj*i£ or within meth- 
ods based on two-component Hamiltonian operators ob- 
tained by other transformation like the Douglas-Kroll- 
Hess transformation^. 

In the proposed combination of spin-orbit resolved 
j-dependent pseudopotentials with the EXX-SCDFT 
method 9 we treat the exchange contribution V x of the 
required exchange-correlation potential V xc exactly and 
neglect the correlation contribution. An inclusion of an 
approximate LDA or GGA correlation potential for the 
component 14,00 would be straightforward. However, ex- 
perience from standard EXX calculations shows that the 
inclusion of an LDA or GGA correlation potential has 
only marginal effect on the bandstructure. 

For the self consistent solid state calculations plane 
wave energy cutoffs of 25 Ry for the orbitals and 14 Ry 
for the KS response function were chosen. The set of used 
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FIG. 1: Bandstructure of germanium including spin-orbit 
coupling 
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FIG. 2: Electron density poo and spin-current pn of germa- 
nium along the bond axis ([111] direction) 

k-points was chosen as a uniform 4x4x4 mesh covering 
the first Brillouin zone. For the lattice constants of silicon 
and germanium the experimantal values of 5.4307A and 
5.6579A, respectively, were used. Throughout our calcu- 
lations, all conduction states were taken into account for 
the construction of the response function and the right 
hand side of the EXX equation^. Pseudopotentials with 
angular momenta / = 0,1,2 and with the j-quantum 
number taking the values j = 1 + 1/2 and j = 1 — 1/2 



were employed for both silicon and germanium. Cutoff 
radii, in atomic units, of ; =0 = 1.8, rf\ =1 — 2.0 and 
r c. \=2 = 2 -° for silicon, as well as r^f =0 = 1.8 rf f =1 = 2.0 
and 7"c f_2 = 2.9 for germanium were used. The germa- 
nium 3<i-electrons were not included in the valence space. 

The resulting EXX bandstructure for germanium is 
shown in Fig. ^ with the values for the indicated spin- 
orbit splittings listed in Table [I] The deviation of the 
calculated SCDFT spin-orbit splittings from the exper- 
imental value is roughly 13.5%. We believe that this 
deviation can be attributed to the fact that for the ger- 
manium calculations the 3<i-states were not included in 
the valence space. This explanation is supported by the 
results for silicon. Here, no d-states are involved and 
the calculated spin-orbit splittings are in almost perfect 
agreement with experiment, see Table [I] 



TABLE I: Spin-orbit splitting for germanium and silicon at 
the T-point 



Germanium 


A(r 7u -r 8 „) 


A(r 6c - r 8c ) 


SCDFT 


258.1 


173.3 


Experiment- 


297 


200 


Silicon 


A(r 25 „) 




SCDFT 


42.5 




Experiment— 


44.1 





Fig. |21 shows the electron density p 00 and the spin- 
current pi2 along the bond axis, i.e., the unit cell's diago- 
nal. The plot for poo shows the usual minimum of the va- 
lence electron density at the positions of the germanium 
nuclei. The currents po^, the spin-densities p^o, and the 
spin-currents p^^ (for p = 1,2,3) turned out to be zero. 
However, the displayed spin-current p\i differs from zero 
and the spin-currents P23 and pj,\ equal p\i and thus also 
are non-zero. This shows that spin-orbit effects indeed 
induce spin- currents. The latter influence the resulting 
bandstructures and therefore should not be neglected in 
a selfconsistent treatment of spin-orbit effects. Further- 
more, the symmetry of the investigated system proposes 
the relation p 12 = p 2 3 = ~Pm = ~p2i = ~P32 = P31 be- 
tween the spin-currents with nonzero value. This relation 
was confirmed by the results of our calculations. 

In summary, we have presented a method to treat non- 
collinear spin, magnetic effects, currents, spin-currents, 
and spin-orbit interactions in solids on an equal footing in 
a self-consistent EXX Kohn-Sham approach and showed 
that spin-orbit coupling induces spin-currents. 



1 R. G. Parr and W. Yang, Density- Functional Theory of 4 see, e.g., A. Gorling, J. Chem. Phys. 123, 062203 (1995); 
Atoms and Molecules (Oxford University Press, 1989). and references therein. 

2 M. Stadele, J. A. Majewski, P. Vogl, and A. Gorling, Phys. 5 S. Sharma et al., cond-mat/0510800 

Rev. Lett. 79, 2089 (1997). 6 G. Vignale and M. Rasolt, Phys. Rev. Lett. 59, 2360 

3 T. Kotani, Phys. Rev. Lett. 74, 2989 (1995). (1987). 



■5 



7 G. Vignale and M. Rasolt, Phys. Rev. B 37, 10685 (1988). 

8 N. Helbig, Diploma Thesis, Universitat Wiirzburg 2001. 

9 S. Rohra and A. Gorling, Phys. Rev. Lett. 97, 013005 
(2006). 

10 E. Engel, A. Hock, and S. Varga, Phys. Rev. B 63, 125121 
(2001). 

11 E. Engel et al., Phys. Rev. B 64, 125111 (2001). 



E. van Lenthe, E. J. Baerends and J. G. Snijders, J. Chem. 

Phys. 99, 4597 (1993). 

B. A. Hess, Phys. Rev. A 33, 3742 (1986). 

O. Madelung, Semiconductors: Data Handbook (Springer, 

Berlin, 2004). 



